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SUMMARY 


The  orbital  parameters  of  Cosmos  220  rocket  (1968-40B)  have 
been  determined  at  74  epochs  from  some  6300  observations.  For  38  of 
the  determinations,  Hewitt  camera  observations  were  available. 

The  orbit  was  determined  from  January  1983  until  July  1985, 
and  during  this  time  the  satellite  passed  through  29:2  resonance. 

The  variations  in  inclination  and  eccentricity  have  been  analysed  to 
determine  six  lumped  geopotential  harmonics  of  order  29.  The  first 
two  have  standard  deviations  equivalent  to  accuracies  of  0.3  and 
0.6  cm  in  geoid  height,  better  than  has  previously  been  determined 
for  order  29  by  resonance  analysis  or  any  other  method.  The  standard 
deviations  of  the  other  four  values  correspond  to  accuracies  between 
1.2  and  3.5  cm  in  geoid  height.  The  results  are  independent  of 
recent  comprehensive  geoid  models  and  are  used  to  test  five  such 
models:  the  comparisons  show  that  the  estimated  errors  given  for  the 
models  are  quite  realistic. 
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1  INTRODUCTION 

The  rocket  which  put  Cosmos  220  into  orbit  on  1968  May  7  itself  also 
entered  orbit  and  was  designated  1968-40B.  Initially  it  had  the  following 
orbital  elements':  inclination  74.05°;  nodal  period  99.17  min;  apogee  and  perigee 
heights,  759  km  and  678  km  respectively;  and  eccentricity  0.006. 

After  some  15  years  in  orbit,  1 968-40B  was  slowly  approaching  the  condition 
of  29:2  resonance,  ie  when  the  satellite  makes  29  revolutions  while  the  Earth 
spins  twice,  relative  to  the  orbital  plane.  As  a  result  1968-4QB  was,  in  1983, 
placed  on  the  priority  list  for  high-priority  observing  and  remained  there  for 
three  years,  so  that  the  orbit  of  the  satellite  could  be  determined  at  frequent 
intervals  while  it  was  passing  through  the  resonance.  In  this  Report  the  orbit 
of  1968-40B  has  been  determined  from  radar  and  optical  observations  between  1983 
January  and  1985  July  using  the  RAE  orbit  determination  program  PROP  in  the 
PR0P6  version2.  The  inclination  and  eccentricity  were  then  analysed  over  this 
time  to  evaluate  29th-order  lumped  harmonics. 

2  THE  OBSERVATIONS,  ORBITS  AND  OBSERVATIONAL  ACCURACY 

2. 1  The  observations 

The  orbit  of  1968-40B  has  been  determined  at  74  epochs,  from  some 
6300  observations,  between  1983  January  5  and  1985  July  27.  The  number  of  obser¬ 
vations  used  in  each  determination  is  given  in  Table  1  and  the  source  of  the 
observations  is  also  indicated. 

There  were  four  groups  of  observations  available.  Those  listed  first  in 
Table  1  were  the  most  accurate,  being  made  by  the  University  of  Aston  s  Hewitt 
cameras  at  the  Royal  Greenwich  Observatory  site  at  Herstmonceux  and  at  the  Siding 
Spring  Observatory  site  in  Australia:  these  observations  usually  have  an  accu¬ 
racy  of  a  few  seconds  of  arc.  The  second  consisted  of  visual  observations  made 
by  volunteer  observers  who  reported  to  the  University  of  Aston,  and  these  obser¬ 
vations  usually  have  accuracies  between  2  and  5  minutes  of  arc  in  good  observing 
conditions.  The  third  group  is  made  up  of  Navspasur  observations  kindly  supplied 
by  the  US  Naval  Research  Laboratory,  with  accuracies  of  about  2  minutes  of  arc. 
The  final  group,  British  radar  observations,  came  from  the  radar  trackers  at 
RAF,  Fylingdales. 

2 . 2  The  orbits 

The  orbits  were  determined  at  approximately  14-day  intervals  with  the 
aid  of  the  RAE  orbit  refinement  program  PROP  in  the  PROP  6  version,  and  the 
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Table  1 


Sources  of  the  observations  used  on  each  orbit 


Orbit 

No. 

Source  of  observation 

— 

Total 

Source  of  observation 

Total 

Hewi 1 1 

camera 

Visual 

US 

Navy 

British 

radar 

Visual 

US 

Navy 

British 

radar 

1 

3 

70 

27 

100 

38 

25H 

10 

25 

21 

81 

2 

1 1 

61 

28 

100 

39 

10H 

7 

18 

12 

47 

3 

1 

66 

33 

100 

40 

13 

41 

36 

90 

4 

12 

61 

26 

99 

41 

13 

50 

31 

94 

5 

12S 

2 

48 

14 

76 

42 

40 

31 

71 

6 

5 

60 

17 

82 

43 

10H 

37 

47 

7 

12H 

3 

12 

56 

83 

44 

10H 

7 

27 

6 

50 

8 

6H 

53 

15 

74 

45 

10H 

21 

21 

17 

69 

9 

10H 

8 

68 

13 

99 

46 

5H 

9 

23 

20 

57 

10 

6H 

64 

6 

76 

47 

1 

36 

33 

70 

1 1 

17H 

7 

47 

29 

100 

48 

2 

43 

38 

83 

12 

6 

47 

31 

84 

49 

8 

40 

34 

82 

13 

1 

55 

44 

100 

50 

5 

54 

38 

97 

14 

2  Oh 

36 

37 

93 

51 

6S 

8 

55 

31 

100 

15 

21H 

3 

31 

26 

81 

52 

24 

49 

27 

100 

16 

1  1H 

3 

26 

23 

63 

53 

6S 

52 

42 

100 

17 

6H 

6 

35 

19 

66 

54 

4S 

65 

26 

95 

18 

5H 

14 

25 

27 

71 

55 

48 

42 

90 

19 

32 

42 

26 

100 

56 

14S 

23 

48 

15 

100 

20 

10 

50 

39 

99 

57 

6S 

4 

48 

40 

98 

21 

56 

36 

92 

58 

4 

45 

51 

100 

22 

3H 

8 

54 

30 

95 

59 

6S 

3 

59 

31 

99 

23 

6H 

2 

36 

20 

64 

60 

21 

52 

27 

100 

24 

6H 

6 

24 

20 

56 

61 

7S 

49 

43 

99 

25 

6H 

31 

18 

55 

62 

9S 

27 

31 

67 

26 

2 

40 

36 

78 

63 

53 

47 

100 

27 

2 

39 

40 

81 

64 

17 

55 

28 

100 

28 

5 

29 

36 

70 

65 

6S 

4 

62 

28 

100 

29 

7 

53 

40 

100 

66 

4S 

1 

39 

36 

80 

30 

9H 

5 

36 

24 

74 

67 

2 

52 

45 

99 

31 

6 

46 

46 

98 

68 

5H 

1 

65 

29 

100 

32 

48 

42 

90 

69 

14 

53 

33 

100 

33 

1 1 

50 

39 

100 

70 

12S 

9 

36 

39 

96 

34 

4 

47 

37 

88 

71 

50 

50 

100 

35 

44 

29 

73 

i  72 

4 

54 

42 

100 

36 

2  OH 

35 

55 

1  73 

8H 

1  1 

35 

42 

96 

I” 

20H 

37 

_ 

10 

67 

74 

4S 

1 

A - 

43 

36 

_ 

84 

i. 

S  =  camera  at  Siding  Spring,  Australia.  H  =  camera  at  Herstmonceux,  England 
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orbital  elements  at  each  epoch  are  listed  in  Table  2,  with  the  standard  deviations 
below  each  value.  The  epoch  for  each  orbit  is  at  00  hours  on  the  day  indicated, 
and  the  PROP  program  fits  the  mean  anomaly  M  by  a  polynomial  of  the  form 

M  =  Mq  +  Mjt  +  M2t2  +  M3t3  +  M^t4  +  M5t5  ,  (1) 

where  t  is  the  time  measured  from  epoch,  and  the  number  of  M  coefficients 
used  depends  on  the  drag.  The  satellite  1968-40B  was  in  a  nearly  circular  orbit 
at  a  height  of  about  650  km,  where  drag  is  slight:  consequently  on  31  orbits  it 
was  only  necessary  to  use  the  coefficients  Mq  and  Mj  ;  on  39  orbits  the 
coefficients  were  needed;  and  on  4  orbits  the  coefficients  Mq  - 

were  required. 

The  value  of  c  ,  the  parameter  indicating  the  measure  of  fit  of  the 
observations  to  the  orbit,  varied  between  0.44  and  1.05,  and  had  an  average 
value  of  0.74.  For  all  74  orbits  the  standard  deviations  in  inclination  were 
between  0.0005°  and  0.0014°,  the  average  being  0.0009°  which  is  equivalent  to 
about  110  m  in  distance.  The  values  of  the  standard  deviations  in  eccentricity 
varied  between  0.000001  and  0.000022,  and  the  average  value  was  0.000011,  corres¬ 
ponding  to  80  m  in  distance.  The  accuracy  for  the  right  ascension  of  the  node 
was  nearly  the  same  as  that  for  the  inclination. 

2. 3  The  accuracy  of  the  observations 

The  331  I  US  Navy  observations  used  in  the  orbit  determination  were  of  the 

usual  consistent  accuracy  of  about  2  minutes  of  arc,  and  those  of  the  British 

radars  were  also  of  their  customary  standard.  The  residuals  of  other  observa- 

3 

tions  have  been  obtained  with  the  aid  of  the  ORES  computer  program  and  sent  to 
observers.  Table  3  lists  the  rms  residuals  of  the  observing  stations  with  five 
or  more  observations  accepted  in  the  final  orbit  determinations. 

The  residuals  of  some  visual  observers  are  a  little  higher  than  usual, 
but  as  these  rms  values  include  all  observations  made  by  each  observer,  with  no 
attempt  to  eliminate  those  made  in  bad  conditions  of  visibility,  this  is  not 
surprising. 

The  Hewitt  cameras  did  very  well:  the  Herstmonceux  camera  contributed 
265  observations  with  rms  residual  of  8  seconds  of  arc,  equivalent  to  about  30  m 
in  position,  and  the  Siding  Spring  camera  contributed  81  observations  with  rms 
residual  of  3  seconds  of  arc,  equivalent  to  about  10  m.  Many  of  the  74  orbits 


Table  2 


Values  of  orbital  parameters  at  74  epochs  with  standard  deviations 


MJD 

Date 

a 

e 

i 

0 

u 

Ml 

Mj 

*3 

C  1 

TT 

45339 

1983  Jan  5 

7067.0870 

0.005515 

74.0474 

28.877 

39.9 

243.858 

5262.0693 

0.0036  1 

0.55  1 

77 

7 

2 

1  1 

12 

1 

1 

3 

1 

45351 

1983  Jan  17 

7066.9854 

0.005118 

74.0476 

5.931 

15.9 

3.338 

5262.1828 

0.0040 

|  0.70 

86 

7 

3 

16 

1  1 

1 

1 

3 

2 

45366 

1983  Feb  1 

7066.8993 

0.004618 

74.0522 

337.247 

343.2 

64.240 

5262.2792 

0.0038 

0.57 

84 

7 

2 

14 

12 

1 

l 

3 

1 

45386 

1983  Feb  21 

7066.7927 

0.004064 

74.0497 

299.006 

293.8 

147.995 

5262.3932 

0.0015 

0.51 

89 

8 

2 

8 

9 

1 

2 

2 

1 

45401 

1983  Mar  8 

7066.7328 

0.004095 

74.0459 

270.318 

258.9 

211.848 

5262.4649 

0.0013 

1.05 

75 

7 

4 

15 

12 

1 

2 

4 

l 

45418 

1983  Mar  25 

7066.6701 

0.004381 

74.0448 

237.805 

211.0 

357.415 

5262.5349 

0.0033 

0.44  1 

65 

9 

2 

10 

8 

1 

1 

2 

1 

45430 

1983  Apr  6 

7066.6081 

0.004786 

74.0477 

214.857 

183.6 

122.248 

5262.6044 

0.0015 

0.47 

73 

7 

45438 

1983  Apr  14 

7066.5794 

0.005067 

74.0497 

199.558 

166.5 

85.849 

5262.6366 

0.0026 

0.46  , 

58 

; 

3 

9 

6 

1 

1 

3 

1 

| 

45455 

1983  May  1 

7066.4942 

0.005597 

74.0491 

167.049 

133.2 

234.604 

5262.7318 

0.0025 

0.00016 

0.56  | 

90 

9 

3 

7 

S 

1 

<  | 

3 

1 

3 

45480 

1983  May  26 

7066.3749 

0.005879 

74.0450 

119.231 

87.0 

350.230 

5262.8649 

0.0035 

0.54 

74 

9 

2 

6 

7 

1 

l 

2 

1 

45496 

1983  Jun  11 

7066.3070 

0.005717 

74.0496 

88.629 

57.6 

282.015 

5262.9410 

0.0022 

0.00029 

0.57 

99 

7 

5 

8 

6 

1 

l 

6 

1 

6 

45508 

1983  Jun  23 

7066.2519 

0.005449 

74.0486 

65.675 

34.9 

51.656 

5263.0025 

0.0047 

0.00057 

0.77 

77 

7 

10 

18 

12 

1 

1 

1  1 

2 

1  1 

45527 

1983  Jul  12 

7066.2061 

0.004838 

74.0459 

29.332 

355.4 

288.039 

5263.0535 

0.0009 

0.69 

98 

; 

2 

12 

6 

1 

1 

2 

1 

45550 

1983  Aug  4 

7066.1584 

0.004214 

74.0482 

345.327 

300. 1 

328.944 

5263. 1069 

0.0008 

0.00028 

0.  74 

91 

1 

9 

8 

1 

1 

1 

2 

8 

45560 

1983  Aug  14 

7066.1272 

0.004092 

74.0531 

326.203 

273.3 

18.466 

5263.1420 

0.89 

81 

•j 

1 

10 

9 

1 

2 

1 

45564 

1983  Aug  18 

7066.  1229 

0.004075 

74.0502 

318.547 

262.0 

182.379 

5263.1468 

1 

0.87 

61 

1 

8 

I  1 

2 

2 

1 

I 

45573 

1983  Aug  27 

7066.0888 

0.004182 

74.0506 

301 .333 

239.5 

11.371 

5263. 1849 

0.66 

63 

7 

10 

8 

1 

2 

8 

i 

45577 

1983  Aug  31 

7066.0803 

0.004267 

74.0492 

293.678 

228.9 

175.454 

5263. 1943 

0.82 

68 

8 

12 

10 

2 

2 

9 

45584 

1983  Sep  7 

7066.0630 

0.004532 

74.0466 

280.289 

211.5 

282.683 

5263.2135 

0.0016 

0.80 

99 

5 

16 

10 

1 

2 

6 

4 

45595 

1983  Sep  18 

7066.0336 

0.004938 

74.0436 

259.240 

186.4 

194.332 

5263.2462 

0.0019 

0.82 

99 

4 

19 

10 

1 

1 

4 

3 

45610 

1983  Oct  3 

7065.9964 

0.005404 

74.0460 

230.543 

154.9 

270.777 

5263.2880 

0.0017 

0.67 

89 

3 

14 

9 

1 

1 

3 

2 

45635 

1983  Oct  28 

7065.9138 

0.005956 

74.0468 

182,716 

107.3 

39.977 

5263.3803 

0.0010 

0.81 

95 

4 

7 

7 

1 

2 

5 

3  : 

45642 

1983  Nov  4 

7065.8915 

0.005982 

74.0479 

169.326 

94.3 

148.519 

5263.4053 

0.0026  1 

0.83 

64 

8 

10 

7 

1 

2 

9 

7 

1 

45646 

1983  Nov  8 

7065.8834 

0.005985 

74.0462 

161 .672 

86.8 

313.466 

5263.4142 

0.0045 

0.69 

55 

7 

9 

6 

1 

1 

7 

7 

45650 

1983  Nov  12 

7065.8670 

0.005974 

74.0464 

154.020 

78.9 

1  18.470 

5263.4326 

0.0050 

0.65 

55 

7 

8 

7 

1 

2 

8 

4 

45662 

1983  Nov  24 

7065.8432 

0.005859 

74.0426 

131.058 

57.2 

253.805 

5263.4590 

0.83 

78 

4 

14 

1  1 

1 

2 

5 

45677 

1983  Dec  9 

7065.8172 

0.005430 

74.0461 

102.358 

28.6 

333.335 

5263.4883 

0.0012 

0.93 

79 

5 

21 

13 

2 

2 

6 

4 

45697 

1983  Dec  29 

7065.7800 

0.004675 

74.0451 

64.089 

345.4 

80. 175 

5263.5298 

1  0.96 

70 

6 

22 

14 

2 

2 

7 

1984  Jan  6 

7065.7664 

0.004469 

74.0462 

48.783 

326.8 

51.101 

5263.5450 

0.0008 

0.83 

97 

4 

13 

10 

1 

1 

4 

2 

45719 

1984  Jan  20 

7065.7425 

0.004236 

74.0441 

21 .992 

291.2 

270.537 

5263.5717 

0.77 

70 

3 

10 

9 

1 

3 

3 

i 

45730 

1984  Jan  31 

7065.7154 

0.004147 

74.0444 

0.942 

261.8 

186.090 

5263.6019 

o.ooio  ; 

0.71 

97 

4 

9 

9 

! 

2 

5 

2 

45748 

1984  Feb  18 

7065.6665 

0.004498 

74.0470 

326.500 

215.7 

212.329 

5263.6567 

0.0013 

0.82 

90 

3 

16 

1  1 

1 

2 

4 

2 

45764 

1984  Mar  5 

7065.6115 

0.005026 

74.0448 

295.884 

179.5 

156.558 

5263.7181 

0.0022 

0.85 

98 

3 

20 

10 

1 

1 

4 

2 

45774 

1984  Mar  15 

7065.5813 

0.005341 

74.0409 

276.748 

158.6 

212.213 

5263.7516 

0.0022 

0.76 

84 

4 

9 

10 

i 

1 

5 

3  ! 

45795 

1984  Apr  5 

7065.4531 

0.005867 

74.0426 

236.553 

117.9 

186.768 

5263.8951 

0.0041 

0.65 

70 

3 

1 1 

10 

1 

1 

4 

2 

45803 

1984  Apr  13 

7065.4144 

!  0.005946 

74.0472 

221.245 

103.3 

160.781 

5263.9385 

0.68 

54 

45812 

1984  Apr  22 

7065.3851 

10.005941 

74.0496 

1204.024 

86.5 

]  356.832 

5263.9714 

0.0019  1 

0.69 

|  66 

<1 

1  6 

7 

I  i 

1 

1 

i 

i  1 
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Table  2  (concluded) 


MJD 

Date 

a 

e 

i 

ft 

u 

u  +  Mq 

Mi 

M2 

*3 

e 

N 

D 

38 

45817 

1984  Apr  27 

7065.3669 

0.005920 

74.0480 

194.453 

77.5 

25.887 

5263.9917 

0.0025 

0.92 

79 

4.0 

1 

6 

8 

1 

2 

1 

l 

39 

45821 

1984  May  1 

7065.353^ 

0.005888 

74.0463 

186.798 

70.0 

193.207 

5264.0068 

0.84 

47 

3.5 

1 

5 

13 

2 

2 

1 

40 

45826 

1 984  May  6 

7065.3397 

0.005820 

74.0426 

177.232 

61.5 

222.425 

5264.0218 

0.0014 

0.90 

89 

5.8  i 

5 

13 

)  2 

2 

2 

5 

3 

[ 

41 

45834 

1984  May  14 

7065.3224 

0.005625 

74.0436 

161.924 

46.2 

197.304 

5264.0412 

0.0014 

0.62 

92 

5.8 

3 

)  1 

8 

1 

1 

3 

2 

42 

45853 

1984  Jun  2 

7065.2512 

0.004985 

74.0409 

125.558 

7.9 

93.605 

5264.1207 

0.0016 

0.69 

70 

5.7 

5 

18 

10 

J 

1 

5 

3 

43 

45861 

1984  Jun  10 

7065.2279 

0.004761 

74.0463 

110.245 

350.0 

69.309 

5264.1470 

0.0022 

0.64 

47 

5.6 

5 

8 

10 

j 

J 

5 

2 

44 

45866 

1984  Jun  15 

7065.2162 

0.004572 

74.0469 

100.680 

338.2 

99.217 

5264.1600 

0.72 

49 

3.7 

2 

8 

1 1 

2 

J 

2 

45 

45870 

1984  Jun  19 

7065.2083 

0.004411 

74.0465 

93.023 

328.7 

267.192 

5264. 1689 

0.80 

66 

3.7 

l 

17 

9 

1 

2 

J 

46 

45878 

1984  Jun  27 

7065. 1916 

0.004262 

74.0477 

77.712 

309.4 

243.255 

5264.1877 

0.73 

53 

3.0 

6 

1 1 

13 

1 

2 

7 

4? 

45897 

1984  Jul  16 

7065. 1647 

0.004077 

74.0422 

41.343 

259.2 

141.867 

5264.2175 

0.87 

70 

5.2 

5 

12 

13 

2 

3 

6 

48 

45906 

1984  Jul  25 

7065.1542 

0.004192 

74.0415 

24.116 

236.0 

340.336 

5264.2291 

0.87 

83 

5.9 

4 

13 

1 1 

1 

2 

5 

49 

45934 

1984  Aug  22 

7065. 1208 

0.005170 

74.0489 

330.524 

170.6 

78.463 

5264.2669 

0.73 

79 

5.3 

4 

19 

10 

1 

1 

4 

i 

50 

45947 

1984  Sep  4 

7065. 1034 

0.005569 

74.0487 

305.644 

144.5 

85.825 

5264.2863 

0.73 

94 

3.8 

3 

14 

9 

1 

1 

4 

51 

45955 

1984  Sep  12 

7065.0897 

0.005768 

74.0445 

290.332 

129.3 

62.806 

5264.3015 

0.0007 

0.59 

98 

3.8 

3 

6 

6 

1 

<1 

4 

2 

52 

45964 

1984  Sep  21 

7065.0783 

0.005916 

74.0426 

273. 106 

112.8 

262.030 

5264.3141 

0.0014 

0.58 

96 

6.2 

3 

8 

7 

I 

1 

3 

2 

53 

45980 

1984  Oct  7 

7065.0527 

0.005940 

74.0446 

242.473 

82.8 

216.544 

5264.3428 

0.0010 

0.60 

98 

6.0 

3 

6 

6 

f 

<1 

3 

2 

54 

45986 

1984  Oct  13 

7065.0438 

0.005878 

74.0482 

230.989 

71.8 

109.605 

5264.3529 

0.0008 

0.60 

69 

5.6 

2 

4 

6 

1 

<1 

3 

2 

55 

46002 

1984  Oct  29 

7065 .0184 

0.005571 

74.0486 

200.369 

41.9 

64.741 

5264.3814 

0.0008 

0.73 

87 

5.6 

4 

14 

9 

1 

1 

5 

3 

56 

46017 

1984  Nov  13 

7065.0017 

0.005091 

74.0460 

171.662 

11.5 

158.028 

5264.3999 

0.59 

97 

5.2 

<  1 

1 

6 

<1 

<1 

1 

57 

46038 

1984  Dec  4 

7064.9680 

0.004338 

74.0415 

131.464 

322.7 

*45.165 

5264.4374 

0.85 

96 

3.8 

4 

17 

9 

1 

1 

4 

58 

46056 

1984  Dec  22 

7064.9546 

0.004092 

74.0504 

97.01 1 

276.6 

186.080 

5264.4527 

0.0003 

0.80 

99 

8.3 

3 

10 

10 

1 

2 

3 

1 

59 

46077 

1985  Jan  12 

7064.9318 

0.004434 

74.0435 

56.811 

222.4 

174.259 

5264.4779 

0.86 

99 

6.5 

4 

6 

10 

i 

1 

5 

60 

46091 

1985  Jan  28 

7064.9176 

0.004829 

74.0463 

30.014 

189.4 

46.650 

5264.4939 

0.0004 

0.57 

96 

6.3 

2 

14 

7 

i 

1 

3 

2 

61 

46104 

1985  Feb  8 

7064.9057 

0.005270 

74.0445 

5.123 

162.0 

56.926 

5264.5072 

0.0008 

0.89 

93 

5.9  1 

? 

5 

5 

1 

1 

2 

3 

i 

62 

46109 

1985  Feb  13 

7064.9007 

0.005423 

74.0473 

355.554 

151.8 

88.619 

5264.5129 

0.86 

67 

4.U  1 

2 

6 

10 

1 

i 

2 

■ 

63 

46121 

1985  Feb  25 

7064.8919 

0.005761 

74.0489 

332.591 

128.4 

236.776 

5264.5229 

0.90 

9' 

4.6  i 

5 

15 

1  1 

] 

2 

6 

64 

46134 

1985  Mar  10 

7064.8821 

0.005964 

74.0468 

307.710 

104.6 

247.432 

5264.5337 

0.89 

97 

5.2 

4 

10 

10 

1 

2 

4 

65 

46144 

1985  Mar  20 

7064.8784 

0.005998 

74.0431 

288.56? 

85.6 

31 1 .077 

5264.4376 

0.63 

90 

5.9 

3 

5 

8 

1 

<1 

3 

i 

66 

46154 

1985  Mar  30 

7064.8742 

0.005865 

74.0417 

269.417 

67.9 

14.763 

5264.5422 

0.93 

78 

6.0 

5 

6 

1 1 

1 

1 

6 

67 

46170 

1985  Apr  15 

7064.8711 

0.005541 

74.0418 

238.788 

37.5 

332.714 

5264.5457 

0.81 

96 

7.5  i 

2 

15 

10 

1 

1 

3 

68 

46182 

1985  Apr  27 

7064.8613 

0.005208 

74.0484 

215.822 

12.8 

J21 .249 

5264.5570 

0.81 

95 

6.0, 

3 

15 

8 

1 

1 

4 

69 

46195 

1985  May  10 

7064.8554 

0.004703 

74.0464 

190. 939 

344.9 

132.293 

5264.5634 

0.56 

95 

i  7.0, 

2 

13 

7 

l 

1 

3 

70 

46205 

1985  May  20 

7064.8464 

0.004390 

74.0462 

171.795 

320. 1 

196. 257 

5264.5735 

0.0012 

0.85 

96 

6.9  | 

2 

6 

8 

1 

1 

2 

1 

71 

46221 

1985  Jun  5 

7064.8382 

0.004106 

74.0433 

141. 164 

279.2 

154.757 

5264.5825 

0.84 

100 

1  7.5, 

3 

9 

10 

1 

3 

4 

72 

46234 

1985  Jun  18 

7064.8289 

0.004175 

74.0450 

1  16.276 

244.3 

166.178 

5264.5930 

0.  68 

98 

7.4  , 

3 

14 

1  1 

1 

2 

4 

73 

46251 

1985  Jul  5 

7064.8232 

0.004576 

74.0489 

83.734 

202.7 

347.410 

5264.5996 

0.88 

96 

6.4  ■ 

3 

18 

8 

1 

l 

4 

74 

46273 

1 985  Jul  27 

7064.8101 

0.005325 

74.0451 

41.614 

154.8 

201.001 

5264.6141 

0.93 

84 

5.3 

5 

5 

12 

1 

1 

6 

KEY 


MJD  modified  Julian  day 

a  semi  major  axis  (km) 

e  eccentricity 

i  inclination  (deg) 

ft  right  ascension  of  ascending  node  (deg) 
m  argument  of  Derieee  (dee) 


Mg  mean  anomaly  at  epoch  (deg) 

M |  mean  motion  n  (deg/day) 

M2, M3  later  coefficients  in  the  polynomial  for  M 

c  measure  of  fit 

N  number  of  observations  accepted  in  each  orbit  determination 

D  time  covered  by  the  observations  (days) 
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Table  3 


Residuals  for  observing  stations  with  five  or  more  observations 
accepted  in  the  final  orbit  determination 


— 

Number  of 

Rms 

| 

residuals 

Station 

observations 

Minutes  of  arc 

Accepted 

Rejected 

RA 

Dec 

Total 

414 

Capetown 

20 

I 

1.7  i 

2.  1 

2.7 

1  184 

Eilenburg 

7 

0 

5.9  ' 

1.9 

6.2 

2265 

Farnham 

24 

3 

5.4  , 

5.0 

7.4 

2392 

Cowbeech 

8 

1 

1.0 

2.4 

2.6 

2414 

Bournemouth 

132 

0 

3.2 

3.5 

4.7 

2418 

Sunningdale 

27 

0 

2.3 

2.2 

3.2 

2420 

Willowbrae 

83 

6 

4.0 

3.4 

5.2 

2430 

Stevenage  4 

8 

0 

1.7 

2.  1 

2.7 

2437 

Warrington 

5 

0 

6.1 

7.4 

9.6 

2539 

Dymchurch 

25 

0 

1.5 

1  .6 

2.2 

2657 

Bridgwater 

30 

0 

1.8 

2.2 

2.8 

2658 

Hillsborough 

8 

0 

3.0 

1.6 

3.4 

2659 

Herstmonceux  3 
(Hewitt  camera) 

265 

2 

0.  10 

0.10 

0.14 

2665 

Cluj  3 

5 

0 

4.3 

5.9 

7.3 

4156 

Apeldoorn 

7 

0 

1.2 

3.8 

4.0 

9652 

Siding  Spring 
(Hewitt  camera) 

81 

_ 

4 

. 

0.03 

0.04 

0.05 

utilize  more  than  one  Hewitt  camera  transit  and  the  PROP  orbital  model  is  not 
accurate  to  better  Liuui  50  m  over  an  interval  of  more  than  a  few  days  because 
lunisolar  perturbations  are  neglected.  So  orbital  error  almost  certainly  contri¬ 
butes  to  the  Hewitt  camera  residuals:  the  cameras  are  more  accurate  than  the 
residuals  indicate.  The  Herstmonceux  camera  suffered  more  from  this  effect, 
because  many  of  the  orbit  determinations  included  multiple  transits:  for  example, 
orbit  38  used  four  Hewitt  camera  plates  at  daily  intervals. 

3  EQUATIONS  FOR  29:2  RESONANCE 

The  theory  for  29:2  resonance  is  detailed  in  Ref  4,  where  all  the  para¬ 
meters  used  are  defined.  The  theoretical  equation  for  the  change  in  inclination 
at  resonance, for  orbits  with  small  eccentricity,  is 
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di 

dt 


29  r 

sin  i  \  a  /  la 


(29  -  2  cos  i)F 

+  16e(29  -  cos  i 

+  I2e(29  -  3  cos  i)F 


30,29, 14) 


(.0,2  _0 , 2 

sin  i>  +  C_.  cos  $ 


29 


29 


i)f29, 29,14^29'  si"(i  -  ">  -  529  C“(t  -  “>j 


_-l,3 


_~'.3  | 

29,29,1 3)C29  Sin(*  +  w)  "  S29  COS(*  + 


q  cos,  .  . 

+  terms  in  e  .  (yv  -  qw) 

sin 


]• 


where 


$  =  2(d)  +  M)  +  29(ft  -  v) 


(2) 

(3) 


is  the  resonance  angle,  u  being  the  argument  of  perigee,  M  the  mean  anomaly, 

fz  the  right  ascension  of  the  node  and  v  the  sidereal  angle.  The  quantities  y 

and  q  are  integers,  but  only  three  terms  are  given  explicitly  in  equation  (2) 

2 

because  it  is  believed  the  others  are  small:  terms  with  q  =  ±2  have  e  as  a 
multiplying  factor  and,  for  I968-40B,  e  is  0.006  at  most  during  the  time  of 
29:2  resonance;  while  the  terms  with  y  =  2  are  associated  with  harmonics  of 
order  58,  which  should  be  much  smaller  than  those  of  order  29.  The  F  in 
equation  (2)  ate  functions  of  inclination  only,  defined  in  Ref  4. 

.  .  .  4 

The  theoretical  equation  for  the  variation  of  e  ,  for  small  e  ,  is 


de 

dt 


,29,I4\ 


'_0 , 2  _0, 2 

29  sin  $  +  C29  cos  0 


) 


.1,1 


-  16F 


+  I2F 


C29  sin!1?  -  w)  -  S29  cos(il>  -  w)| 


29 


,29, 13{ 


-1,3 


-1,3 


C29  sin ($  +  w)  -  S29  cos ($  +  w)! 


terms  in  £e^  '|q  -  (y  +  9)e "  |  gi^^  ~  9W) 


]] 


(4) 


The  main  terms  in  equation  (4)  are  expected  to  be  those  with  (y,q)  =  (1,1)  and 
(1,-1),  but  the  (y,q)  =  (1,0)  term  is  included  so  that  the  analysis  conforms  with 
that  for  inclination. 

_0, 2 

The  quantities  such  as  S29  in  equations  (2)  and  (4)  are  lumped  geo¬ 
potential  harmonics  of  order  29,  which  can  be  expressed  as  linear  sums  of 

individual  coefficients  C„  and  S„  , 

im  m 


10 


r 


m 


and 


(5) 


and  these  individual  coefficients  can  be  evaluated  when  enough  lumped  harmonics 
trom  satellites  at  different  inclinations  are  available.  The  summation  for  l 
in  equation  (5)  is  in  steps  of  2,  beginning  at  the  lowest  l  ,  which  is  either  29 
or  30  here  (see  Ref  4),  Also  k  =  2y  -  q  for  29:2  resonance. 


The  orders  of  magnitude  of  the  terms  in  the  equations  for  di/dt  and 

de/dt  can  be  estimated,  as  the  C 

52  -q,k 

H  ,  so  the  value  of  C  (or 

— 5  2  2 )  i  ® 

x  10  J/l  )  J  ,  the  Q  coefficients  being  obtained  for  specific  values  of 
(y,q)  using  the  RAE  computer  program  PROF. 
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(or  S.  )  are  expected  to  be  of  order 
_q,k 

S  )  can  be  taken  to  be  of  order 
m 


Estimating  the  orders  of  magnitude  of  the  terms  in  equation  (2)  for 
inclination  suggests  that,  if  just  the  (y,q)  =  (1,0)  term  is  used,  the  likely 
errors  are  21  from  the  neglect  of  the  (2,0)  term  and  8%  from  the  neglect  of  the 

(1,±1)  terms.  In  equation  (4)  the  (y,q)  =  (1,0)  term  only  contributes  about 

.  2 
0.5%.  As  mentioned  before,  the  q  =  2  terms  contain  a  factor  e  and  therefore 

they  will  not  be  needed. 


4  ANALYSIS  OF  INCLINATION  AT  29:2  RESONANCE 

Cosmos  220  rocket  passed  through  exact  29:2  resonance  on  1984  Aug  15  and 
its  orbit  has  been  determined  either  side  of  resonance,  with  the  rate  of  change 
of  resonance  angle,  $  ,  increasing  from  -4.3  to  +0.7  deg/day.  The  variation  of 
both  $  and  the  resonance  angle  $  ,  given  by  equation  (3),  are  shown  in  Fig  1. 
The  decrease  of  4>  became  very  slow  in  the  last  18  months  due  to  the  reduction 
in  solar  activity. 


Before  the  changes  in  inclination  due  to  resonance  can  be  analysed,  all 
other  known  perturbations  must  be  removed.  The  74  values  of  inclination  in 
Table  2  were  therefore  cleared  of  lunisolar  and  zonal  harmonic  perturbations,  by 
using  the  computer  program  PROD"*  with  I -day  integration  steps,  and  the  perturba¬ 
tion  due  to  the  2  tesseral  harmonic,  which  is  recorded  on  each  PROP  run,  was 
also  removed  from  each  value  of  inclination.  Perturbations  due  to  earth  and 
ocean  tides  should  not  exceed  50  m  and  need  not  be  considered,  since  the  values 
of  inclination  have  an  average  accuracy  of  110  m,  the  most  accurate  being  0.0005° 
which  is  approximately  equivalent  to  65  m. 
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The  values  of  inclination,  cleared  of  the  perturbations  mentioned,  and  with 
the  standard  deviations  quoted  in  Table  2,  were  fitted  with  equation  (2)  in  inte¬ 
grated  form  using  the  RAE  THROE^  computer  program.  This  program  removes  the  per¬ 
turbations  due  to  atmospheric  rotation  and  lunisolar  precession  of  the  Earth's 
axis  before  fitting  the  remaining  perturbation  due  to  resonance.  The  density 

scale  height,  H  ,  was  taken  as  85  km,  appropriate  to  a  height  of  683  km,  0.4H 

7  8 

above  perigee  ,  and  the  atmospheric  rotation  rate,  A  ,  was  taken  as  0.9  rev/day. 

The  values  of  were  altered  to  mean  values,  ,  by  the  technique  described 

in  Ref  9. 

The  first  fitting  of  the  values  of  inclination  by  THROE,  with  (y,q)  =  0,0), 
gave  c  =  2.29  ,  where  e  is  the  measure-of-f it  parameter.  On  inspection  of  the 
PROP  runs  for  some  of  the  ill-fitting  values  of  inclination,  it  was  found  that 
the  runs  included  quite  a  few  low-elevation  radar  observations,  and  relaxing  the 
standard  deviations  of  the  range  values  had  a  beneficial  effect  on  the  values  of 
inclination.  After  these  relaxations  the  value  of  e.  on  the  next  THROE  run  was 
I.C8)  and  the  values  of  the  lumped  harmonics  were  as  follows: 

0  2  0  2 

109C2g  =  -2.6  ±0.5  ,  1  °9S2g  =  8.9  ±1.0  .  (6) 


For  this  fitting,  one  value  of  inclination  had  its  standard  deviation  increased 
by  a  factor  of  1.5  to  ensure  that  all  weighted  residuals  were  less  than  2c  . 

The  addition  of  the  (y,q)  =  (1,±1)  terms  did  not  improve  the  standard 
-  -  0  2 

deviation  of  (C,S)09  and  the  additional  harmonics  were  not  well  determined.  The 
values  obtained  from  the  THROE  run  were: 


_0, 2 

C29  =  -  2.2  ±  0.5 


15  ±  13 


-1,3 

C29  =  -  94  ±  36 


(7) 


with  e  =  0.927  .  The  values  of  (C,S)  ’  and  (C,S)?  ’  were  of  the  right  order 
-5  2  ^ 

according  to  the  10  /&  rule,  but  they  were  poorly  determined.  The  results  are 
given  because  the  run  was  needed  for  the  simultaneous  fitting  of  inclination  and 
eccentricity.  The  addition  of  the  (y,q)  =  (2,0)  terms  also  proved  unhelpful. 


12 


The  values  of  inclination,  cleared  of  all  known  perturbations  except  those 
due  to  resonance,  are  plotted  in  Fig  2.  The  theoretical  curve  derived  from  the 
THROE  fitting  with  (y,q)  =  (1,0),  that  gave  the  values  in  equation  (6),  is  also 
shown,  as  a  full  line. 

5  ANALYSIS  OF  ECCENTRICITY  AT  29:2  RESONANCE 

Before  the  changes  in  eccentricity  due  to  resonance  can  be  analysed,  the 
values  of  eccentricity  in  Table  2  have  first  to  be  cleared  of  perturbations  due 
to  zonal  harmonics,  atmospheric  drag  and  lunisolar  perturbations. 

The  removal  of  zonal  harmonic  and  lunisolar  perturbations  was  performed  by 
using  the  PROD  computer  program,  and  the  variation  due  to  drag  in  an  atmosphere 
with  day-to-night  variation  was  calculated  using  the  theory  given  in  the  Appendix 
of  Ref  10  with  scale  height  H  taken  as  106  km,  appropriate  to  a  mean  height  of 
734  km.  This  theory  assumes  that  the  atmospheric  density  depends  on  the  geocen¬ 
tric  angular  distance  from  the  point  of  maximum  density,  which  has  been  taken 
as  14  h  local  time*^. 

In  Fig  3a  the  values  of  eccentricity  from  Table  2  are  given  as  circles,  and 
these  values  after  removal  of  zonal  harmonic  and  lunisolar  perturbations  are 
indicated  by  triangles.  The  correction  AeQ  for  drag  in  an  atmosphere  with  day- 
to-night  variation  is  shown  in  Fig  3b.  The  values  of  eccentricity  cleared  of 
these  perturbations  are  now  ready  to  be  fitted  with  equation  (4)  in  integrated 
form,  using  THROE. 

After  the  values  of  had  been  altered  to  mean  values,  ,  as  for  the 

inclination,  the  first  THROE  fitting  with  (y,q)  =  ( 1 , ± 1 )  was  performed  and  the 
value  of  e  obtained  was  3.91.  However,  two  of  the  eccentricity  values,  had 
standard  deviations  less  than  0.000005;  so  both  were  increased  to  that  value,  in 
case  there  was  a  significant  effect  from  the  neglect  of  earth  and  ocean  tides. 

Six  values  of  eccentricity  had  to  have  their  standard  deviations  increased  by  a 
factor  of  2,  and  two  by  a  factor  of  4,  to  keep  all  the  weighted  residuals  less 
than  2e  .  The  THROE  run  with  these  adjustments  gave  e  =  2.53  . 

In  some  previous  fittings  of  values  of  eccentricity  with  THROE,  it  was 
found  that  adjustment  of  the  odd  zonal  harmonics  was  needed.  The  same  method  has 
been  used  here  as  in  Ref  12,  ie  to  add  an  increment  AJ^  to  the  value  used 

in  the  PROD  model.  The  lowest  value  of  e  ,  2.497,  was  obtained  with 
AJ3  =  -0.015  x  10"6  . 
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This  final  fitting  of  the  values  with  (y,q)  =  0,-0  and  with 
AJ^  =  -0.015  x  10  6  ,  gave  values  of  the  lumped  harmonics  as  follows: 

9-0 '  9-01 

10  C29  =  53.4  ±  3.8  10  s29  =  _10-3  1  5,6 

q_-03  q  -1  ,3 

10  C2g  =  18.9  ±  6.9  10  S  =  ~3*9  1  7*8 

The  values  of  eccentricity,  cleared  of  all  known  perturbations  larger  than  the 
standard  deviations  of  the  values  (except  those  due  to  resonance),  are  plotted 
in  Fig  4.  The  fitted  theoretical  curve  is  that  of  the  THROE  run  which  produced 
the  values  in  equation  (8),  and  is  shown  as  a  full  line. 

6  INCLINATION  AND  ECCENTRICITY  FITTED  SIMULTANEOUSLY 

The  values  of  inclination  and  eccentricity  fitted  separately  by  THROE  can 
be  fitted  simultaneously  using  the  RAE  computer  program  SIMRES  developed  by 
Dr  R.H.  Gooding.  This  program  combines  the  results  from  a  number  of  THROE  runs 
(each  with  the  same  set  of  (y,q)  terms),  and  produces  a  single  set  of  coefficients 
to  fit  the  data.  For  this  SIMRES  fitting,  the  results  of  THROE  runs  with 
(y,q)  =  (1, 0)(1,1)  and  (1,-1)  were  used.  The  SIMRES  program  allows  a  choice  of 
weighting,  so  that  the  contributing  THROE  runs  can  be  given  more  or  less  weight 
according  to  their  accuracy  of  fit,  which  is  indicated  by  the  value  of  e  . 

The  THROE  fitting  of  inclination  with  (y,q)  =  (1,0) (1,1)  and  (1,-1)  gave 
e  =  0.927  ,  and  for  eccentricity  the  value  of  c  was  2.100,  when  fitted  with  the 
same  terms.  For  the  SIMRES  fitting,  therefore,  the  weighting  of  e  was  down¬ 
graded  by  a  factor  equal  to  the  ratio  of  the  values  of  e  on  the  THROE  fittings, 
namely  2.265  (=2.100/0.927).  The  values  of  the  lumped  harmonics  given  by  this 
SIMRES  fitting  are: 

q_°  >  2  q_0,2 

1°  C29  "  ~2*7  1  °-5  10  S29  =  8’2  1  K° 

9J.l  9_1,1 

,OC29  =  51.2  ±3.5  >0S29  -  -7.5  ±5.3 

q  -1,3  Q-"1’3 

10C29  =  16.1  ±6.5  10S29  =  -2.3  ±7.3 

The  SIMRES  fittings  of  inclination  and  eccentricity  are  shown  pictorially  in 
Figs  2  and  4,  as  broken  lines. 
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The  lumped  harmonics  obtained  from  the  fitting  of  inclination  and  eccen¬ 
tricity  separately,  given  in  equations  (6)  and  (8),  are  very  similar  to  those 
obtained  from  the  simultaneous  fitting,  given  in  equation  (9):  all  like  coef¬ 
ficients  are  within  one  standard  deviation  of  each  other.  So  either  set  of 
values  could  be  used  in  future  determinations  of  the  individual  29th-order 
harmonic  coefficients  in  the  geopotential.  The  SIMRES  values  are  probably  better, 
for  two  reasons.  First,  they  take  into  account  the  (y,q)  =  ( 1 , ± 1 )  terms  for 
inclination,  and  these  may  be  significant  because,  when  estimating  the  orders  of 
magnitude  of  the  terms  (see  section  3),  it  was  found  that  neglect  of  the  (1,±1) 
terms  in  fitting  the  inclination  could  give  an  error  of  8%.  Second,  the  SIMRES 
curve  does  seem  to  provide  the  better  fit  in  Fig  2,  especially  between  MJD  45650 
and  45900,  and  e  is  slightly  lower  (1.05  as  against  1.08). 

7  EQUATIONS  FOR  INDIVIDUAL  COEFFICIENTS 


The  lumped  harmonics  in  equation  (9)  can  be  expressed  as  linear  sums  of 

29th-order  individual  coefficients:  see  equation  (5).  The  Q  coefficients  in 
.  4 

these  equations  depend  on  the  ratios  of  the  eccentricity  functions  »  a°d 

can  be  evaluated  using  the  RAE  computer  program  PROF.  However,  a  correction 

factor  has  to  be  applied  to  the  values  of  Q  ,  as  it  was  assumed  in  the  PROF 
2 

program  that  the  e  terms  in  the  functions  could  be  neglected  (see 

Ref  4).  This  leads  to  large  errors  in  Q  if  e  is  large;  but  here,  with  e 
approximately  0.005,  the  correction  is  almost  negligible,  being  less  than  1%. 

The  correction  has,  nevertheless,  been  made.  The  values  of  the  lumped  harmonics 
from  THROE  and  SIMRES  also  need  a  very  small  correction,  which  has  been  made  in 
equations  (6)  to  (9). 


The  equations  given  below  for  the  individual  coefficients  have  been  termi¬ 
nated  when  the  expected  contribution  from  the  coefficients  permanently  falls  to 
less  than  5%  of  the  largest  contribution. 

The  resulting  six  equations  are  as  follows: 


^29  =  C30,29  +  0*044S2,29"  0,356S4,29'0'349S6,29"0'175^38,29  °*000C40,29 


+  0-l01C42,29  +  0'122C44,29  • 

_°,2 

S^q  :  the  equation  is  the  same  as  (10)  with  S  instead  of  C  , 


(10) 


ri  n 
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J.l 

C29  =  C29,29  "  1  -  576C31  29  "  0>726C33,29  +  °‘244S5,29  +  °*  635S?  ,29  +  °‘  524^39,29 

+  °-19854.,29  -0-099543,29  -  0-249S45,29  -°-246SA7,29  ’  °-148S9,29‘ 


S„_  :  the  equation  is  the  same  as  (12)  with  S  instead  of  C  , 


C29  ‘  C29,29-0'590C3I,29  -  0-759C33,29-0-464C35,29-  0'087C37,29  *  0',83C39,29 


*°-290C4,,29  *  °'258C43,29  *  °'I48C45,29  * 


:  the  equation  is  the  same  as  (14)  with  S  instead  of  C  .  (15) 

8  LUMPED  HARMONICS  FROM  1968-40B  COMPARED  WITH  THOSE  FROM  COMPREHENSIVE  CF.OID 
MODELS 


It  is  interesting  to  compare  the  values  of  lumped  harmonics  obtained  here 

with  those  from  comprehensive  gravity-field  models.  Five  models  have  been  used 

1  3 

for  comparison:  GEM  10B  and  GEM  10C  ,  the  model  produced  in  ’981  by 

R.H.  Rapp14,  GRIM  3— L 1  *  ^  and  GEM-T1  18  .  The  values  of  29th-order  lumped  har¬ 
monics  given  by  substituting  the  values  of  the  individual  coefficients  from  the 
models  into  equations  (10)  to  (15)  appear  in  Table  4. 


The  Goddard  Earth  Model  GEM  I  OB  extends  to  order  and  degree  36,  and  GEM  10C 
consists  of  the  GEM  10B  solution  up  to  degree  36,  together  with  some  31000  coef¬ 
ficients  of  order  and  degree  up  to  180,  derived  from  analysis  of  altimeter 
measurements  over  the  oceans.  Equations  (10)  to  (15)  require  coefficients  above 
36,  some  as  high  as  degree  49,  but  the  lumped  harmonics  in  Table  4  for  GEM  10B 

are  truncated  at  degree  36.  The  expected  accuracy  of  the  GEM  10B  coefficients 

1 7  -9 

has  been  assumed  to  be  3  «  10  .  The  GEM  10C  accuracy  (above  degree  36)  is 

1 8  -9 

certainly  poorer  and  is  rather  arbitrarily  taken  as  5  x  10  .  The  standard 

deviations  for  the  GEM  10B  and  GEM  10C  values  in  Table  4  have  been  assessed 
using  these  accuracies  for  the  individual  coefficients. 


The  comprehensive  geopotential  model  produced  by  R.H.  Rapp  at  the  Ohio 
State  University  in  1981  gives  the  individual  coefficients  to  order  and  degree 
180,  and  is  derived  from  Seasat  altimeter  data,  terrestrial  gravity  measurements 
and  other  data.  An  accuracy  estimate  for  each  coefficient  is  also  given,  so  the 
standard  deviations  for  the  lumped  harmonics  in  Table  4  have  been  assessed  using 
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these  accuracies.  Finally,  the  lumped  harmonics  obtained  by  using  the  individual 
coefficients  from  GRTM3-L1  and  GEM-T1  are  given  in  Table  4.  These  models  only 
extend  to  degree  and  order  36,  so  the  equations  (10)  to  (15)  have  to  be  truncated 
at  this  value,  as  with  GEM10B.  The  standard  deviations  are  calculated  from  the 
accuracies  given  for  the  individual  coefficients  in  both  cases. 


Table  4 


Values  of  lumped  harmonics  from  1968-40B  and  comprehensive  geoid  models 


9-°’ 2 
10C 

29 

q_0,2 

10  S 

29 

qj*1 

>°c29 

qj.l 

.0  s29 

o.-l, 3 
10C29 

Q_~  1 >  3 

10S29 

1 968-40B 

2.7  ±  0.5 

8.2  ±  1.0 

51.2  ±  3.5 

-7.5  ±  5.3 

16.1  ±  6.5 

-2.3  ±  7.3 

GEM  10C 

-0.5  ±  3 

5.2  ±  3 

25.9  ±  8 

-12.5  ±  8 

15.3  ±  5 

-5.2  i  5 

GEM  10B 

1 .3  ±  3 

4.9  ±  3 

24.7  ±  6 

-8.9  ±  6 

14.4  ±  4 

-5.4  ±  4 

Rapp  (1981) 

-0.5  ±  1.9 

7.1  ±  2.5 

32.0  ±  4.2 

-1.2  ±  6.5 

13.3  ±  3.1 

-11.7  ±  4.3 

GRIM3-L! 

0.8  ±  3.3 

5.8  ±  3.3 

34.5  ±  6.0 

-10.3  +  6.0 

24.6  ±  4.3 

-13.6  ±  4.3 

GEM-T1 

6.5  ±  6.2 

2.7  ±  6.2 

34.7  ±10.4 

11.6  ±  10.4 

24.2  ±  8.8 

8.0  +  8.8 

On  examination  of  Table  4  it  is  seen  that  all  the  lumped  harmonics  from 

—  1  9  1 

the  models,  except  C 29  ,  agree  well  with  those  from  1968-40B,  within  0.6  times 

—0 , 2 

the  sum  of  their  standard  deviations  on  average;  only  the  Rapp  value  for  Con 

-1,1  29 

and  the  GEM-T1  value  for  differ  from  the  I968-40B  values  by  more  than  the 

29  _1,1 

sum  of  their  standard  deviations.  The  values  for  do  not  agree  so  well 

with  that  from  1968-40B;  on  average  the  difference  is  2.1  times  the  sum  of  the 

-1.1 

standard  deviations.  The  value  of  from  1968-40B  is  large  and  well-defined, 

originating  from  the  large  and  well-defined  increase  in  eccentricity  between 
MJD  45400  and  45700  in  Fig  4. 

The  standard  deviations  of  the  values  from  1968-40B  are  much  better  than 
those  of  the  models  for  the  first  pair  of  coefficients,  somewhat  better  for  the 
second  pair,  and  rather  worse  for  the  third  pair. 

The  good  agreement  between  the  lumped  harmonics  in  Table  4  is  very  satis- 
fying,  as  it  shows  that  the  29th-order  coefficients  in  the  models  must  be  of  the 
right  magnitude  and  that  their  estimated  errors  are  quite  realistic.  A  previous 
analysis  of  the  29:2  resonance  of  1 96 7  —  1 04B  with  less  accurate  orbits  led  to 
similar  though  less  definite  conclusions.  Hopefully  in  the  future  more  satellites 
can  be  analysed  at  other  inclinations  as  they  pass  through  29:2  resonance,  and 
then  the  individual  coefficients  can  be  evaluated  from  satellite  orbit  analysis. 
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In  conclusion,  it  is  worth  emphasizing  that  the  results  obtained  here  from 
1 968-40B  are  completely  independent  of  the  models:  the  gravity  field  used  in 

PROP  is  from  a  much  earlier  model,  of  zonal  harmonics  and  J  only;  and 

z ,  / 

I 968-40B  is  not  among  the  satellites  utilized  in  the  models  of  Table  4. 

9  APPROXIMATE  ACCURACY  IN  GEO ID  HEIGHT 

Equations  (10)  to  (15)  are  useful  in  allowing  an  approximate  assessment  of 
the  accuracy  of  the  lumped  harmonics  evaluated  for  1968-40B,  in  terms  of  an 
accuracy  in  geoid  height.  The  accuracy  o  in  geoid  height  may  be  estimated 
approximately  as  Ro/Q*  ,  where  a  is  the  error  in  the  lumped  harmonic  and 


The  values  of  for  each  lumped  harmonic  are  as  follows: 

_0, 2  _0 , 2  J,  1  _1,1  _-l  ,3  1,3 

^29  ^29  ^29  ^29  ^29  S29 

a  (cm)  0.3  0.6  1.2  1.8  3.1  3.5 

8 

The  values  for  (q,k)  =  (1,1)  are  more  accurate  than  those  for  (q,k)  =  (-1,3), 

probably  because  the  orbit  determination  covered  the  (y,q)  =  (1,1)  resonance 

•  • 

(<t>  =  w)  ,  but  did  not  continue  long  enough  to  cover  fully  the  (y,q)  =  (1,-1) 

•  • 

resonance  (4>  =  -u) .  The  accuracy  of  the  first  pair,  (q,k)  =  (0,2),  is  much 

better  than  has  previously  been  obtained  for  order  29  from  resonance  analysis,  or 
any  other  method  (as  Table  4  shows). 

10  CONCLUSIONS 

The  orbit  of  1968-40B  has  been  determined  at  74  epochs  from  some  6300 
observations,  between  1983  January  and  1985  July  while  the  satellite  was  passing 
through  the  condition  of  29:2  resonance.  The  average  accuracy  of  the  inclination 
and  eccentricity  for  all  74  epochs  was  equivalent  to  110  m  and  80  m  in  distance 
respectively. 

The  variations  in  inclination  and  eccentricity  have  been  analysed,  and  six 
29th-order  lumped  harmonics  have  been  evaluated:  the  recommended  values  are  those 
given  in  equation  (9).  The  first  two  have  standard  deviations  equivalent  to 
accuracies  of  0.3  and  0.6  cm  in  geoid  height,  considerably  better  than  has 
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previously  been  obtained  for  order  29  from  resonance  analysis  or  any  other  method. 
The  standard  deviations  of  the  other  four  values  correspond  to  accuracies  between 
1.2  and  3.5  cm  in  geoid  height.  The  results  provide  an  independent  test  of  the 
29th-order  harmonics  from  recent  comprehensive  gravity  models:  comparisons  with 
five  such  models  (Table  4)  show  good  agreement  and  suggest  that  the  estimated 
errors  of  the  models  are  quite  realistic,  for  order  29. 
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